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ABSTRACT 



This annual report details progress in basic research in thermoacoustic 
heat transport made during the period October 1, 1991 through September 30, 
1992. Our major research efforts were concentrated in two areas during FY 1992: 
1) a standing wave analysis of the low amplitude performance of thermoacoustic 
prime movers both below and above the onset of self-oscillation and 2) a study of 
finite amplitude standing waves in both harmonic and anharmonic tubes. In 
addition to these areas, FY 1992 research efforts also included the investigation of 
heat driven refrigerators and preliminary work on investigations of 
thermoacoustic engines with laser Doppler anemometry. This technical report 
summarizes the important results of these projects. A publications, patents, 
presentations, and honors report is also included. 
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SUMMARY OF PROGRESS 



Our major research efforts in thermoacoustic heat transport were 
concentrated in two areas during FY 1992: 1) a standing wave analysis of the low 
amplitude performance of thermoacoustic prime movers both below and above 
the onset of self-oscillation and 2) a study of finite amplitude standing waves in 
both harmonic and anharmonic tubes. In addition to these areas, FY 1992 
research efforts also included the investigation of heat driven refrigerators and 
preliminary work on investigations of thermoacoustic engines with laser Doppler 
anemometry. 



1) Standing Wave Analysis of the Low Amplitude Performance of Prime 
Movers 

A. Below Onset of Self-Oscillation 

The purpose of this project is to develop a method of analyzing the 
performance of prime movers and comparing the predictions to experimental 
results. This analysis, motivated by Swift, 1 assumes that the acoustic pressure 
and velocity can be described in terms of a standing wave throughout the prime 
mover and is not limited by either the small boundary layer or short stack 
approximations. It falls somewhere between the two methods of analysis 
discussed in Ref. 1. It is more general than the short stack, boundary layer 
approach but not as general as the fully numerical approach. It maintains the 
conceptual simplicity of the former and allows one to appreciate the need for the 
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later. The discussion that follows summarizes a complete discussion of the 
analysis found in Ref. 2. 

The results of the analysis are shown for the first longitudinal mode of a 
helium filled prime mover in Fig. 1. The figure, a graph of the reciprocal of the 
quality factor of the prime mover Q versus the temperature difference applied 
across the prime mover stack, corresponds to a mean gas pressure of 170 kPa. 
The solid line indicates the results of our standing wave analysis. For 
comparison purposes, the results of the porous medium analysis (discussed in 
Ref. 3 and the result of our FY 1991 research) is included, the dashed lines. It is 
evident that the standing wave analysis provides a better fit to the 170 kPa data 
than does the porous medium approach. However, the two techniques do more- 
or-less equally well at higher mean gas pressures. There is a tendency for the 
standing wave technique to under predict the attenuation, which is proportional 
to 1/Q, at zero temperature difference. 

The data presented in this figure are a composite of individual data sets 
recorded at different times with (slightly but inevitably) different mean 
temperatures and mean gas pressures. The spread in the data is indicative of the 
uncertainty. Sources of experimental errors are uncertainties in mean 
temperature, temperature difference, mean gas pressure, the positions of the 
prime mover elements, and the geometry (plate spacing and perimeter) of the 
heat exchangers and prime mover stack. The net effect of these error sources is 
illustrated using the results presented in Fig. 1. The predictions were 
recalculated including each individual error source and the deviations between 
the new and the original calculations found. These individual deviations were 
added together to arrive at the total deviation. In Fig. 2, the solid line represents 
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the original calculation as presented in Fig. 1. The dashed (dash-dot) line 
represents the original calculation plus (minus) the total deviation due to errors 
in mean temperature, temperature difference, mean pressure, and positions and 
perimeters of the elements. Most of the data points fall within these limits, 
though there are some exceptions, especially near onset. These errors introduce 
an uncertainty in the temperature difference required for onset of approximately 
20 °C. 



The calculations discussed up to this point have been somewhat involved. 
The heat conduction equation is solved numerically, taking into account the 
temperature dependence of the thermal conductivity of the gas and stainless steel 
stack plates. An example of the distribution of the temperature and temperature 
gradient is shown in Fig. 3 for a temperature difference of 350 K. It is seen that 
the temperature (solid line) is almost linearly distributed along the stack. The 
temperature gradient (dashed line) varies by approximately 20% over the length 
of the stack. This temperature distribution is used to find the spatial dependence 
of the temperature dependent variables, to allow numerical integration along the 
stack. 



A natural and important question arises. How much can these procedures 
be relaxed and still obtain reasonable results? Providing, at least, some of the 
answers to this question is an important result of this work. Neglecting the 
thermal conductivity of the gas has negligible effect on the results of the 
calculations. Assuming that the temperature gradient is constant along the stack 
(and simply equal to the temperature difference divided by the length of the 
stack) has only a slight effect. It pushes the onset temperature higher than that 
predicted by the full calculation by about 5 °C at 170 kPa. The difference 
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diminishes as the mean pressure increases, virtually disappearing at 500 kPa. All 
of these simplifications are probably tolerable unless one requires very precise 
results. 

A much more drastic simplification is to assume a constant temperature 
gradient and use the average stack temperature to calculate the thermophysical 
properties of the gas along the entire length of the stack. This assumption greatly 
simplifies the analysis, removing the need for numerical integration. At 170 kPa, 
this assumption reduces the predicted onset temperature by approximately 
10 °C. However, at 500 kPa the situation is somewhat worse. In this case, the 
predicted onset temperature is reduced by approximately 25 °C. Therefore, 
while this assumption certainly has computational advantages, it does not 
produce very accurate results near onset. However, for temperature differences 
sufficiently below onset, the differences between predictions based on this 
assumption and the full calculation are likely tolerable. 

To test our ability to predict the Qs of higher modes, we measured the Q 
as a function of temperature difference for the second and third modes of the 
prime mover. The results for a mean gas pressure of 170 kPa are shown in Fig. 4. 
The solid lines show the predictions of the standing wave analysis. As is the case 
with the fundamental mode, the analysis tends to under predict 1/Q. However, 
the overall trends match well. The dashed lines are the results of correcting the 
predicted values by adding to them the differences between the predicted and 
measured results at zero temperature difference. The corrected prediction agrees 
quite well with the data, especially for temperature differences well below onset. 
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B. Above Onset - During the Initial Growth of Oscillations 



The comparisons presented above were made for a prime mover 
operating below onset. Below onset the temperature distribution in the stack is 
governed by thermal conduction through the plates. Under these circumstances 
the analysis gives good agreement with the data. However, above onset, the 
steady state acoustic pressure amplitudes are typically 1 - 10% of mean pressure. 
At these high amplitudes the temperature distribution in the stack is significantly 
influenced by the acoustic heat transport. It remains to be seen how well the 
present analysis predicts the large amplitude performance of a prime mover 
operating above onset. Of course, given the actual temperature distribution in 
the stack, there is no reason to expect that the standing wave analysis would 
break down. However, if one is going to go the effort of accurately calculating 
the temperature distribution under large amplitude conditions, then one might 
as well as do a fully numerical solution to the problem. However, it is difficult to 
gain intuition from a fully numerical approach unless one has a fair amount of 
experience with either thermoacoustics or numerical calculations. It is suggested 
that, because it is conceptually relatively simple, the standing wave analysis is a 
good alternative. 

Low amplitude conditions do exist in prime movers above onset - during 
the initial build up of the oscillations. By measuring the rise time of the buildup, 
Q can be determined and plotted on the same graph with the below onset Qs 
determined from frequency response measurements. Such data are shown in 
Fig. 5, for a helium filled prime mover at a mean gas pressure of 376 kPa. In 
contrast to the earlier figures, this time we plot Q, not 1/Q, vs. AT. In this 
representation, onset is clearly defined by the rapid divergence of Q. The solid 
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line is the result of the standing wave analysis. The agreement is quite good over 
a 400 K span of AT. (It is worth pointing out that there are no adjustable 
parameters in the model.) Therefore, the conclusion drawn earlier, that the 
standing wave analysis does a good job of predicting low amplitude prime 
mover performance, is valid both above and below onset. 



C. Below Onset Revisited 

Practically all thermoacoustics research has concentrated on the large 
amplitude regime, the emphasis being on industrial and commercial applications 
of heat engines. However, Fig. 5 points out another potential use for 
thermoacoustic prime movers - but below onset. Notice that the quality factor of 
the resonator (after all prime movers are resonators) can be adjusted over a range 
spanning at least one order of magnitude. This ability might be useful in 
applications requiring high Q acoustic resonators, such as in photoacoustic 
spectroscopy. 

Consider Fig. 5 once again. What happens to Q if AT is negative? It goes 
down. This point is illustrated in Fig. 6, the only difference between it and Fig. 5 
being that the calculations have been extended to negative values of AT. For this 
prime mover, Q can be reduced by a factor of approximately 2 from its value at 
AT = 0 by applying a AT of -200 K, achievable with liquid nitrogen. This 
reduction may not be all that impressive, but on the other hand this prime mover 
is not optimized for the purpose of reducing Q! However, if the stack was 
optimized (for instance, by using a larger stack surface area to compensate for the 
limited available temperature difference), there is no reason why the Q could not 
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